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This paper deals with the problems of consistence and strong consistence 
of the maximum likelihood estimators of the mean and variance of the drift 
fractional Brownian motions observed at discrete time instants. A central 
limit theorem for these estimators is also obtained by using the Malliavin 
calculus. 

1. Introduction. Long memory processes have been widely applied to various 
fields, such as finance, hydrology, network traffic analysis and so on. Fractional 
Brownian motions are one special class of long memory processes when the Hurst 
parameter H > 1/2. The stochastic calculus for these processes has now been 
well-established (see [2]). When a long memory model is used to describe some 
phenomena, it is important to identify the parameters in the model. In this paper, 
we shall consider the following simple model 

Y t = f it + aBt I , t>0, (1.1) 

where /i and a are constants to be estimated from discrete observations of the 
process Y. Our method works for fractional Brownian motions of all parameters. 
So in this paper we assume that {Bt*,t > 0) is a fractional Brownian motion 
of Hurst parameter H G (0, 1). But we do not discuss the case H = 1/2, the 
standard Brownian motion case since it is known. This means, (B^,t > 0) is a 
mean Gaussian process with the following covariance structure: 

E (Bf Bf )) = \ {t 2H + s 2H -\t- s\ 2H ) . 

We assume that the process is observed at discrete time instants (ii , t 2 , • • • , t/v)- 
To simplify notation we assume tk = kh, k = 1, 2, • • • , N for some fixed length 
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h > 0. Thus the observation vector is Y = (Y tl ,Y t2 , ■ ■ ■ ,Y tN )'. We will ob- 
tain the maximum likelihood estimators fiN and a 2 N of // and a 2 respectively and 
study their asymptotic behaviors. In particular, the almost sure convergence and 
the central limit type theorem. 

The first reason we chose to study dl.ll ) is because it is simple and we can obtain 
explicit estimators. The second reason is that it is also widely applied in various 
fields. The logarithm of a widely used geometric fractional Brownian motion, 
which is popular in finance, is of the form (11.11) . This paper is also complementary 
to the work [6], where the parameter estimation problem (with continuous time 
observation) for fractional Ornstein-Uhlenbeck processes is studied. 

The parameter estimation problem for long memory processes have been well- 
studied (see [1], [3], [4], [9], [11]). Although most work requires the process to be 
stationary, we may still adapt their idea to analyze above model dl.lb . But we shall 
use the method of [6] which seems to be the simplest one to us. This method is 
based on a result of ( [8]) and uses the idea of Malliavin calculus. 

We introduce notation 

Y = /it + aBf, (1.2) 

where and for the rest of the paper t = (h, 2h, ■ ■ ■ , Nh)' and Bf = (B^ , ■■■ , B$ h ) 
The joint probability density function of Y is 

h(Y) = (27ra 2 )-^|rH|-5exp(-^(Y- /U t) / r w 1 (Y-M) 
where 

r* = [[Cov[Bg,Bf h] l^^ N = (f" + f" - \i - j? H \ m .., N • 

The maximum likelihood estimators of [i and a 2 from the observation Y are given 
by 

A = — K i~, (1-3) 

t'r^t 

~ 2 1 (Y / r^ 1 Y)(tT^ 1 t) - (t'P^Y) 2 

N t'T'H 

In Section 2, we shall show that fx and a 2 converge to fi and a 2 both in mean square 
and almost surely. In Section 3, we prove central limit type theorem. In Section 4, 
we give some simulation to demonstrate our estimators /} and a 2 . 

2. Consistence. In this section we will consider the L 2 consistency and the 
strong consistency of both MLE fj, and a 2 . 

Now, let us first consider the L 2 consistency of dl.31 ). 
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Theorem 2.1 The estimator fi (defined by M.3\) ) of \i is unbiased and it converges 
in probability to \i as N — > oo. 

PROOF. Substituting Y by fit + cxBf in (fI31) . we have 

p, = H + a ^-r-^. (2.1) 

tT^t 



Thus E [fi] = /i and hence jj, is unbiased. On the other hand, we have 
Var[/i] = a 2 E 



2r tT^Bf (Bf)T^ti _ 2 tT^r^t _ a 2 



(tT^t)2 



tT^t) 2 tT^t 



Denote 



M = (mij)ij = i t _ tN , where m ij = ]^(i 2H +j 2H -\i-j\ 2H ) 1 
and denote by m^~- the entry of the inverse matrix M~ l of M. Then we may write 

rr 2 rr 2 rr 2 h- 2H ~ 2 

Var[ul = = h~ 2H hr 2 



We shall use the following inequality (with x = N = (1, 2, • • • , N)) 

II x II 2 
x'M" 1 ^ > " 112 , 

"•max 

where A max is the largest eigenvalue of the matrix M. Thus we have 



Var[/t] < cx 2 /i 



|| N » 2 



12 

Since ||N|| 2 = l 2 + 2 2 + . . . + n 2 = " (n+1) 6 (2n+1) we know that || N || 2 « N 3 . 
On the other hand we have by the Gerschgorin Circle Theorem (see [7], Theorem 
8.1.3) 

N 

Xmax < | mij \< CN 2H+1 , 

'-' j= i 

where C a positive constant whose value may be different in different occurrences. 
Consequently, we have 

Var[/x] < Ca 2 h- 2H - 2 N- 3 N 2H+1 = CN 2H ~ 2 . 

which converges to zero as N — > oo. ■ 
Next we study the estimator a 2 defined by (11.41) . 
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Theorem 2.2 We have 

E (^ ) = ^_V aM y^N^Q (22) 

PROOF. By replacing Y with /it + aBp in dT~4b , we have 



n 1 t*r~H 



Thus 



l j N i\ t j h t tfrzh 

ri 



<r 2 ,, r tT^BpffBfMr; 1 *, iV-1 , 
= iv (iv fr^'t ) = — " ■ <23) 

To compute the variance of a 2 we also need to compute E[(a 2 ) 2 ]: 



E\(a 2 ) 2 ] 



a 



1 



N 2 



.//Vr-inH ( t/r H lB f) 2 \2 



((Bf)T^B 



^[((Bf )T^Bf ) 2 ] - 2£[(Bf )T^B^ 



=^ (^[((Bf )l#Bf ) 2 ] - 2E[(Bf)T H 1 Bf (t ^ 1 ^ )2 ] + 3) . (2.4) 

Denoted = r^ 1/2 Bf . ThenE (XX 1 ) = E(T~ 1/2 Bf (Bf )T^ 1/2 ) = J. There- 
fore, X is a standard Gaussian vector of dimension N. For any A small enough and 
e £ E let us compute the following . 

E[exp(A(Bf)'r^Bf + e tT^Bf)] = £[ex P (A|X| 2 + etT^X)] 



(27T) 

A standard technique of completing the squares yields 

£[exp(A(Bf )T^Bf + £ tT^Bf )] = (l-2A)~f exp { |^f^} =: /(A,e) • 
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We are only interested in the coefficient of A 2 and Xe 2 in the above expression 
/(A,e). We have 



1 + ^— (l + 2A + ---) + - 



f(X,s) = (1 + NX + N(N + 2)A 2 + • 

= 1 + NX + N(N + 2)A 2 + • ■ ■ + (JV + 2)Xe 2 t'F- 1 t + ■■■ . 
Comparing the coefficients of A 2 and Ae 2 we have 

E [(Bf)T^Bf) 2 ] = N(N + 2), (2.5) 

^(Bfyr^Bftt'r-iBf) 2 ] = (iv + 2)(t'r- 1 t). 

Hence, we have 

t'r^Bf) 2 (N + 2)(t'T H h) 



EUBPyTT^B?- — / ' ] = K — ' ~> =N + 2. (2.6) 



Using M3\) . UM . (1231 ) and UM, we obtain 



Var[<x 2 ] = E[(a 2 ) 2 ) - {E[a 2 ]) 2 

; [N(N + 2) - 2(N + 2) + 3 - (N - 1 



7V2 1 

1 <r 4 , (2.7) 



2(iV-l) , 



AT2 

which is convergent to 0. Thus we prove the theorem. ■ 
Now we can show the strong consistence of the MLE fi and a 2 as N — > oo. 

Theorem 2.3 77ze estimators fi and a 2 defined by 61. 3D an<i rti.4D , respectively, are 
strongly consistent, that is, 

fx — > ll a.s. as N — > oo (2.8) 
a 2 — > a 2 a.s. as N — > oo (2.9) 

PROOF. Let's prove the convergence for /} first. We will use a Borel-Cantelli 
lemma. To this end, we will show that 

J2P(\V-V\ > ^) <oo (2.10) 



N>1 



for some e > 0. 
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Take < e < 1 — H. Then from the Chebyshev's inequality and the Nelson's 
hypercontractivity inequality [5], we have 

P(|A-Ml>^) = N^E{\^-^)<C p N^{E{\(i-^)Y 12 
< c' p a p h~ ( - 2H+2 ^ p N 2pe+ ( 2H ~^ P . 

For sufficiently large p, we have 2pe + (2H — 2)p < — 1. Thus ( 12.101 ) is proved, 
which implies (12.81) by Borel-Cantelli lemma. 

In the same way, we can show (12.9I ). ■ 



3. Asymptotic. Now we are interested in the central limiting type theorem for 
the estimators fi and a 2 . First from ( 12. 11 ), it is easy to see that 

y / t*r^t(/l - n) ^ JV(0, cr 2 ) as N tends to infinity 

We want to study a 2 
Theorem 3.1 We have 

(a 2 -a 2 ) ^ M(0,a 4 ) as N -> oo . (3.1) 

PROOF. To simplify notation we assume H > 1/2. The case H < 1/2 is 
similar. We define 

r l^(n 2 n 2 \ ° 2 fm^Vr-iR^ (tT^Bf) 2 [W , 

Gn = ] J^ -a) = - / =[(B t )T H B t - ttr _ h ]- ] J-a . 

From (12.7I ). it is obvious that E[G 2 N ] converges to cr 4 . Thus from Theorem 4 of [8] 

to show (13.11) . it suffices to show that || DGn ||^ L2 ^\ q 
First, using the definition of Malliavin calculus, we obtain 

DsG N = ^ - 2 [^(B? )'^Bf - ^ lB ^ W ] , 
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where L> s (Bf )' = (l[ 0| h](a), l[o,2h](«)> ■ • • > 1 [o,iV/ 1 ](s))- Therefore, we have 



II D s Gn \\h 
2a 4 f [' 
'— aH Uo 



\2H-2 



t*I#t 



ill 



D u (Bf)'T- . 
2d 4 



t'l^t 



duds 



N 



4a H 




\u 



\2H~2\ 



D s {Bf)'Y- H l Bf ■ D w (Bf)T^B 



o jo 



t*r-H 

f (tT^.tTg^.tT^ B Bf ^ 



(t*r^t 



4r/i(i) 



=2cr 4 [A 



2A {2) + A 



(3)] 



Since both D S (B^ )T^ 1 Bf and A,(Bf )T^ 1 Bf are Gaussian random variables 
we can write 



2 

! AT2 



/ EiD^T^B* ■ D r (B? )T- 1 Bf ] 

7[0,T] 4 

■ E[D u (Bt)'T^l Bt ■ D v {B?)'v£B?]\s - u 
A / [D.iBfYrgDrB* • D U (B?)'T^D V B? 

iV 7[0,T1 4 



2iJ_2 |r - v\ 2H - 2 dsdrdudv 



s — u r — t; ~ dsdrdudv . 



Let T^ 1 = (r. • 1 )ij=i,...,Ar , = (1^)^=1,...^ and <5 Zfc be the Kronecker 



symbol. We shall use fj^ h | s 



\2H-2 



dsdu = Tut and Y,j=i r i/ r i'j = S U'- 
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Then we have 



E(\A$ -EA$\ 2 ) 

2 f _ x 

= ^ / ^o.wjry ipjfciCn • i[o,i'ft](«)r i /-/i[oj'ft](«) 

iV J[0,T] 4 J 

•a/f|s — u| ~ ajf|r — v| ~ dsdrdudv 

2 N N 

ij=l i',j'=l 

2 N 2 
= — V tf 2 = — 

jV 2 ^ ij ' N ' 
i,i'=i 



which converges to as N — ► oo. 

1(2) r/1 (2)|2N 



Now we deal with . 



s(i4^-E4 2 f) 

2a 



N^kn^ s{Bt) » Bt H 1 " t'r^t 



t*r-H 

• \s — v\ ~ \u — r\ ~ dvdsdudr 

" ^ 7[o,t]4 t'r^t 

. ^(Bf^t-t^Z^Bf _ 2H _ 2 _ r fH-2 dvdsdudr 

_2_ r Ejj=i Eiy^i jjM^jgV^ • ^Hy 1 ^'/^") 
t*i#t 

• Off|s — u| ~ a,H\u — r\ ~ dvdsdudr 
2 

= ]^2 ^jhr^fi'hr-HhT-^k'hr^Tj^ , 

where the summation is over 1 < i,j,i',j',k,l,k',l' < N. Sum first over 1 < 
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hj' < N and then over 1 < I', k < N, we have 

2 Ef,i, k ',i'=xjhlhk'hl'hr^} 2 



iV 2 



A^2 • 



<l -> oo. 

As for we have 

;4 3) | 2 

t'r^Bf • t'r^DyBf t'r^^Bf • t'r^ivBf 



which converges to as iV — > oo 

I (3) 



s(l4 3) -£;4 3) | 2 ) 



2a| 

■ AT2 



[0,T] 4 



(t*r^t)2 



(t*r^t)2 



• (^(t'r^Bf) 2 ) 2 ^- 

■ AT2 



|2#-2|„/ „./|2,H"-2, 



[0,T] 4 



(t*r^t) 2 

|2_H"-2|„/ „,/|2.H-2, 



(t*r^t 



(t'r-H) 2 ^- 



2o| 

■ N 2 

2 " 

^ - 



t'r /i - 1 z) s B^t / r i:i - 1 L' tt B^ 



[0,T] 2 



t*r^t 



Is — ul 2 ^ 2 dsdu 



[0,T] 2 



t*r-H 



a#|s — u\ 2H 2 dsdu 



2 rEfi=i E^'=i ihT-^i'hTJ,), 



N 2 



t*r^t 



1 J -n 

1 «' 



2 _ 2 

~ iV 2 " 



which converges to as N — > oo. 

By triangular inequality, we have that 



2 \2 



£7(|| #Gyv \\ z n -E || DGtv |&) 

-_ E (A® + 4 2) + 4 3) - s(4 1) + 4 2) + 4 3) )) 2 



<9[ j b(4 ) - ^(4 ij )) 2 + £(4 j - ^(4 j )) 2 + £(4 j - £(4 3j )) 



I (2) 



(2)>x2 



,(3) 



l(3)xx21 



This completes the proof of the theorem. 
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4. Simulation. This section contains numerical simulations of the estimators 
obtained in this paper. The fractional Brownian motions are simulated by the Pax- 
son's method [10]. 

TABLE 1 



The means and standard deviations of estimators (fi=0.7880, a 2 =0.8116) 





H= 


=0.25 


H= 


:0.45 


H= 


=0.55 


H= 


:0.75 




A* 


a 2 


A* 


a 2 


A* 


a 2 


A* 


a 2 


Mean 


0.7862 


0.8152 


0.7884 


0.8153 


0.7911 


0.8126 


0.7678 


0.7910 


Std.dev. 


0.0116 


0.0830 


0.0112 


0.0937 


0.0514 


0.0692 


0.0974 


0.0736 



TABLE 2 

The means and standard deviations of estimators (fi=1.5880, a 2 =1.8116) 





H= 


:0.25 


H= 


:0.45 


H= 


=0.55 


H= 


:0.75 




A* 


a 2 




a 2 


A* 


a 2 


A* 


a 2 


Mean 


1.5863 


1.8694 


1.5882 


1.8719 


1.5961 


1.8647 


1.5925 


1.7864 


Std.dev. 


0.0148 


0.1724 


0.0456 


0.1786 


0.0710 


0.1567 


0.1879 


0.1644 



TABLE 3 

The means and standard deviations of estimators (fi=3.5880, a 2 =5.8116) 





H= 


:0.25 


H= 


:0.45 


H= 


=0.55 


H= 


:0.75 






a 2 




a 2 




a 2 


A* 


a 2 


Mean 
Std.dev. 


3.5861 
0.0314 


5.8133 
0.1648 


3.5810 
0.0792 


5.8192 
0.1737 


3.5837 
0.0905 


5.8229 
0.1031 


3.5834 
0.0526 


5.8346 
0.1026 



From these numerical computations, we see the estimators are excellent both 
fori? > 1/2 and H < 1/2. 
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